Systematic Study of Triaxial Deformation 
in the Relativistic Mean Field Theory 

S. Sugimoto'^''''[], K. Sumiyoshi"'|], D. Hirata'^''='|, B.V. Carlson'^'g, 
1. Tanihata"'[|, and H. Toki^'^'Q 

"The Institute of Physical and Chemical Research (RIKEN), 
Hirosawa, Wako, Saitama 351-0198, Japan 
^Research Center for Nuclear Physics(RCNP), Osaka University, 
Mihogaoka, Ibaraki, Osaka 567-0047, Japan 
Japan Synchrotron Radiation Research Institute (JASRI) - SPring-8, 

Sayo, Hyogo, 679-5198, Japan 
^Departamento de Ffsica, Instituto Tecnologico da Aeronautica - CTA, 
Sao Jose dos Campos, Sao Paulo, Brazil 



^ e-mail: satoru@post man . riken . go . j p 
^e-mail: sumi@postman.riken.go.jp 
^ e-maik daisy @sp8sun . springS . or . j p 
^e-maik brett@fis.ita.cta.br 
^ e-maik t anihat a@ rikaxp . riken . go .j p 
^e-mail: toki@rcnp.osaka-u.ac.jp 



1 



Abstract 

We use the relativistic mean field (RMF) theory to systematically study the change 
of deformation of even-even nuclei in the proton-rich Xe region. We investigate the ap- 
pearance of triaxial deformation in 25 nuclei in the region covering Z = 50 — 58 and 
A'^ = 64 — 72 by performing constrained, triaxially symmetric RMF calculations of their 
energy surfaces. We include pairing correlations using the BCS formalism. We find that 
the Sn isotopes are spherical and the Te isotopes are very gamma unstable with shallow 
minima around 7 = 60°. Adding more protons, the Xe , Ba and Ce isotopes have pro- 
late deformations with their sizes increasing with proton number. The neutron number 
dependence is found to be small. We compare the calculated results with the available 
experimented data on the binding energy and the radii. 
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1 Introduction 



The recent progress of radioactive nuclear beam facilities has provided us with marvelous find- 
ings in nuclear physics. Exotic structures such as neutron halos [|I| and neutron skins 0, Q have 
been found in experimental studies of light unstable nuclei in the neutron-rich region. Much 
new information on the shapes and structures of nuclei far from stability is being revealed by 
the systematic measurement of radii and moments of unstable nuclei @, ^. Planned facilities 
in the world will access a large number of unstable nuclei in the whole region of the nuclear 
chart and enable us to explore where and how exotic phenomena of nuclear structure appear 
in the region far from the stability line |^. One of the great interests is to know where the 
deformation of unstable nuclei appears and how the shape of these nuclei changes along the 
isotopic and isotonic chains. 

At the same time, the relativistic many-body framework has been extensively applied to 
study nuclei and nuclear matter p, p. This has been motivated by the recent success of 
the relativistic Brueckner-Hartree-Fock (RBHF) theory, in which the strong density- dependent 
repulsion arises automatically from the relativistic many-body treatment, in reproducing the 
saturation property of nuclear matter [10, 11, 12 1. Among other properties, the relativistic 
mean field (RMF) theory, which is the phenomenological framework of the RBHF theory, has 
been shown to be excellent at describing the properties of unstable nuclei as well as stable ones 



1131, 14, 15, 16, 17|. The RMF theory has also been successfully applied to the study of the 



deformation of nuclei, as well as other properties of the stable and unstable nuclei ilq, 20 



Furthermore, the RMF theory has been used to calculate the equation of state(EOS) of nuclear 
matter in the wide density and temperature regions tabulated for the application to supernova 
simulations . Recently, a systematic study of all even-even nuclei up to the drip lines in the 
nuclear chart has been performed in the RMF theory with axial deformation The ground 
state properties of about 2000 even-even nuclei from Z = 8 to Z = 120 have been studied 
and all possible deformations of each nuclide have been surveyed using a constrained, axially 
deformed RMF model. Through the systematic analysis of the ground state deformations thus 
found, the pattern of the appearance of prolate and oblate deformations has been obtained. 
In the same study, it was also found that the coexistence of prolate and oblate shapes with 
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similar binding energies occurs in many nuclei in the nuclear chart. This coexistence suggests 
the possible appearance of deformation beyond the axial kind, such as triaxial or even higher 
order multipole deformations. 

The appearance of triaxial deformation in this context has been studied in the case of the 
neutron- rich Sulfur isotopes . In the axial RMF calculation for neutron-rich Sulfur isotopes, 
the energy curves as a function of the (3 deformation have two minima, at both prolate and 
oblate deformations, with energies very close to each other. Judging solely from the energy 
curve, one cannot conclude which ground state deformation is realized or whether yet another 
type of deformation appears. The RMF calculations with triaxial deformation of the same 
isotopes have been performed to clarify this point and a smooth shape transition from prolate 
to oblate shapes through triaxial shape has been found along the Sulfur isotopic chain pSf . This 
example motivates us to study further the appearance of triaxial deformation in other regions of 
the nuclear chart. It is interesting to explore where triaxial deformation appears in the nuclear 
chart, especially in relation with the behavior of the appearance of axial deformation. 

In the present study, we have chosen to explore the proton-rich Xe region for the appearance 
of triaxial deformation. We have made a systematic study of 25 even-even nuclei covering 
Z = 50 — 58 and = 64 — 72, using the RMF theory with triaxial deformation, in order to 
clarify how their shapes change as a function of N and Z in this region. We have calculated the 
energy surface of those nuclei as a function of the deformation parameters, (3 and 7, to explore 
the ground state deformation. The previous study of 54Xe, 55CS and seBa isotopes using the 



RMF theory with axial deformation [0] was successful in reproducing the general features of 
the ground state properties. However, disagreement with the measured isotope shift for the 
proton-rich region, which might be due to triaxial deformation, was observed. In the systematic 
RMF calculation with axial deformation which we will discuss in Sect. 3, the shape change 
from oblate to prolate shape occurs as Z increases, in a region in which the two shapes coexist. 
Thus, the axially symmetric RMF calculations strongly suggest that this region could contain 
triaxial deformed nuclei. 

This region has been discussed as a possible region for triaxial deformation in studies using 



conventional frameworks |25, 26, 27, 2q|. There have also been experimental efforts to measure 
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excitation energies in order to study the collective nature of the nuclei in this region |2^, pO| . 
Studies of triaxial deformation within the mean field approach have also been performed in 



other regions of the nuclear chart [31, 32, 33, 34]. 

In Sect. 2, we describe the framework of the RMF theory with deformation. We discuss the 
behavior of the shape within the RMF theory under the assumption of axial symmetry in Sect. 
3. We present the results of the calculations in the RMF theory with triaxial deformation in 
Sect. 4. Results of the calculations are discussed in Sect. 5. We summarize the paper in Sect. 
6. 



2 Relativistic mean field theory 



We briefly describe the framework of the RMF theory and the procedure of the calculation. All 



details can be found in p, |14|, In the RMF theory, the system of nucleons is described by 
fields of mesons and nucleons under the mean field approximation. We start with the effective 
lagrangian, which is relativistically covariant, composed of meson and nucleon fields. We adopt 
a lagrangian with non-linear a and uj terms ||35|| , 



RMF 
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where the notation follows the standard one. On top of the Walecka a - uj model with pho- 
tons and isovector-vector p mesons, non-linear a meson terms are introduced to reproduce the 
properties of nuclei quantitatively and give a reasonable value for the incompressibility P^]. 
The inclusion of the non-linear term of uj meson |]3^ is motivated by the recent success of 
the relativistic Brueckner-Hartree-Fock theory [0, Deriving the Euler-Lagrange equations 
from the lagrangian under the mean field approximation, we obtain the Dirac equation for the 
nucleons and Klein-Gordon equations for the mesons. The self-consistent Dirac equation and 
Klein-Gordon equations are solved by expanding the fields in terms of harmonic-oscillator wave 



functions [|l4|, 



The RMF model contains the meson masses, the meson-nucleon couphng constants and the 
meson self-couphng constants as free parameters. We adopt the parameter set TMA, which 
was determined by fitting the experimental data of masses and charge radii of nuclei in a wide 
mass range [^2|, ^ The parameters are listed in Table 1. We remark that this parameter 
set has a mass dependence so as to reproduce nuclear properties quantitatively from the light 
mass region to the superheavy region. With the TMA parameter set, the symmetry energy 
is 30.68 MeV and the incompressibility is 318 MeV. Note that the bulk properties of nuclear 
matter at saturation with the parameter set TMA is calculated for uniform matter in the limit 
of infinite mass number. 

The TMA parameter set has been used for the systematic study of all even-even nuclei 
up to the drip lines in the nuclear chart within the RMF framework under the assumption of 
axial symmetry p2|. It has been shown that the overall agreement of the calculated results 
using TMA with the experimental data of masses and charge radii is excellent and is found 
to be much better than the results of spherical RMF calculations with TMA. In the present 
study, we extend the RMF calculation with the TMA parameter set from the axially symmetric 
to the triaxially symmetric one, in order to explore the appearance of triaxial deformation in 
the region in which axial deformations with similar binding energies coexist. We discuss the 
correspondence between the axial and triaxial RMF calculations in the subsequent section. 

In order to take into account triaxial deformation, the fields are expanded in terms of 
the eigenf unctions of a triaxially deformed harmonic oscillator potential p3|. We perform 
calculations that constrain the quadrupole moments of the nucleon distribution, in order to 
survey the coexistence of multiple shapes and to identify the ground state deformation. We use 
a quadratic constraint to calculate a complete map of the energy surface as a function of the 
deformation |D|. We take a basis of up to = 12 major shells of the harmonic oscillator 
wave functions. This is normally enough for the constrained calculations in this mass range 
p3| . We have performed calculations with = 14 major shells in the case of no-pairing and 
found that the energy surface does not change significantly for moderate deformation. We have 
also performed non-constrained calculations with pairing up to A^ = 16 and the convergence 
was generally good with A^ = 12. 
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We take the pairing window as given by Gambhir et al. [T^ as £ — A < 2(41A -'^/^)MeV. 
As for the pairing correlations, we perform the RMF calculations with triaxial symmetry using 



a BCS formalism Since we calculate the energy surface in the full /3 — 7 range, we take 
the pairing interaction strength for a given nucleus as G = 23/^4 [Me V] for both protons and 
neutrons ^ . 



Although the BCS-type treatment has often been used in the RMF calculations, it would be 
preferable to incorporate the pairing correlations in the relativistic many body framework in a 
consistent manner. A study with the relativistic Hartree-Bogoliubov theory has been performed 
under the assumption of spherical symmetry [Q. A systematic study of nuclear deformation 



within such a relativistic many-body framework is currently being made |44 



3 Numerical results 

Before we present the results of the RMF calculation with triaxial symmetry, we discuss the 



corresponding results of the RMF calculation with axial symmetry |2^. We examine the 
behavior of the calculated ground state properties of even-even nuclei in the proton-rich Xe 
region. We note here that the calculated properties of deformations as well as masses and radii 
in this region agree very well with the experimental data. 

In the region of Z = 50 — 58 and = 64 — 72, roughly speaking, the deformation changes 
according to the proton number, with a few exceptions. The Sn isotopes {Z = 50) are dominated 
by a spherical shape and the Te isotopes {Z = 52) have, for the most part, an oblate shape. 
The Xe, Ba and Ce isotopes (Z = 54 — 58) all have a prolate shape. We show the energy 
surfaces for Te and Xe isotopes as functions of (3 deformation (Fig. 1). The general trend 
shows a transition in shape from spherical, through oblate, to prolate as the proton number 
increases. 

Furthermore, there is always more than one energy minimum and shape coexistence is 
generally observed in this region of nuclides, as can be seen in Fig. 1. For each nucleus 
there is a corresponding second minimum which has the deformation parameter, j3, of opposite 
sign to that of the absolute minimum. The energy difference between the absolute minimum 
and the second minimum is generally small in this region. 



We next present the RMF calculations with triaxial deformation for 25 nuclei covering the 
range Z = 50 — 58 and = 64 — 72. We present the energy surfaces of these nuclei in the 
f3 and 7 deformation parameter space, calculated using constraints on the two deformation 
parameters. Figure 2 displays the energy surfaces of the Sn, Te, Xe, Ba, Ce (Z = 50 — 58) 
isotopes with N = 64 — 72, arranged in the form of the nuclear chart. The spacing of contours 
is 1 MeV in total energy in all figures. The energy minimum is marked by the black region, in 
which the energy difference is less than 1 MeV from its absolute minimum energy. 

As for the Sn isotopes, the minima appear consistently at the spherical shape [jS = 0). 
The Te isotopes show a 7-soft character, having similar energies along the 7 direction. Most 
of the Xe, Ba and Ce isotopes have minima at prolate deformation. In some cases (^^"^Ba, 
^^^Ba, ^^^Ce, ^^^Ce and ^^''Ce), the region of the shallow minimum extends to quite large 7 
deformation. These results indicate that the triaxial shape is not stable in those nuclei, but 
the energy surfaces are very 7 soft. 



4 Discussion 

We discuss here the influence of the pairing correlations and the effective interaction on the 
triaxial RMF calculations. Since the consistent calculation of pairing correlations with deforma- 
tion within the relativistic many-body framework is still under development, we have performed 
the RMF calculations with pairing correlations in the BCS formalism, as a first study of the 
appearance of axial and triaxial deformation. We see here the effect of pairing correlations on 
the magnitude of the binding energy and the deformation by comparing the results of calcula- 
tions with and without pairing. We show in Fig. 3 the energy surfaces in the /3 — 7 plane of 
^^°Te, ^^^Xe and ^^'^Ba without and with pairing correlations. The qualitative behavior of the 
two cases is very similar. Generally speaking, the magnitude of the (3 deformation is reduced by 
the pairing correlations. The energy minimum seen at finite 7 deformation in ^^^Ba is washed 
out by the inclusion of the pairing correlations. 

In Fig. 4 we show the proton deformation parameters (3 extracted from the RMF calculation 
with axial deformation. In the same figure we also show the experimental data obtained from 



the B(E2) valuespH]. In this figure we see that our results compare well with experimental 



data. 

We also calculate the energy surface of ^^^Ba with alternative NLl parameter set [1^ of the 
RMF theory in order to test the dependence of triaxial deformation on the choice of the RMF 
parameter. In Fig. 5, we compare the energy surfaces obtained by using the TMA and NLl 
parameter sets. A well-distinguished minimum is seen at prolate deformation for the case of 
the NLl parameter set. This feature is slightly stronger than in the TMA case. We mention 
that triaxial deformation is not found in calculations within the Skyrme-Hartree-Fock(SHF) 
theory PT|. It would be interesting to compare the energy surfaces of RMF and SHF theories 
in a wide mass range. 

In Fig. 6 we show the total binding energy of the nuclei studied. The theoretical values 
are taken from the axially symmetric calculations, since no distinguished triaxial shapes were 
found in the triaxial calculations. We see slightly over binding for the Sn isotopes, which may 
be due to the use of the constant pairing correlations and will be studied further. We show in 
Fig. 7 the charge radius as a function of the neutron number. The general tendency is found 
to be quite satisfactory. 



5 Summary 

We have studied systematically the triaxial deformation of 25 even-even nuclei in the proton-rich 
Xe region. We have calculated their ground state structures in the RMF theory with triaxial 
deformation and with pairing correlations and obtained their energy surfaces in the plane of the 
deformation parameters, (3 and 7, by constraining the quadrupole moments. We have explored 
the appearance of triaxial deformation in the region covering Z = 50 — 58 and = 64 — 72, 
by looking for the minima of the derived energy surfaces in the triaxial deformation parameter 
space. Through comparisons with the results obtained in the RMF calculations with axial 
symmetry, we have discussed the correspondence between the coexistence of axial shapes and 
the appearance of triaxial shapes. 

We have found no distinguished energy minima at triaxial deformations. However, the en- 
ergy surfaces are often very 7 soft. This feature is caused by the pairing correlations since, 
when we remove the pairing correlations in the RMF calculations, we find well distinguished 
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triaxial deformation in this mass region. We have compared the energy surfaces of two param- 
eter sets, the TMA and NLl ones. The TMA parameter set provides more softness in the 7 
direction than the NLl one provides. Comparisons of the binding energies and deformations 
with experimental data are, in general, quite satisfactory. 

We note here that we have not worked out the angular momentum and particle number 
projection, which have already been developed in the non-relativistic approach. The restora- 
tions of these symmetries may change somewhat the results on the deformation, as has been 
discussed in the non-relativistic description of deformed nuclei. The relativistic approach to 
triaxial deformation is not yet at such a level of systematic study and these refinements have 
not been considered here. This is certainly a direction for future work. 
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Table 1 



niN [MeV] 


938.900 


rria [MeV] 


519.151 


rriu, [MeV] 


781.950 


rrip [MeV] 


768.100 


g-j 


10.055 + 3.050/A°-^ 




12.842 + 3.191/^°-^ 




3.800 + 4.644/^0-4 


92 


-0.328 - 27.879/^0-4 




38.862 - 181.191/_4"-4 


C3 


151.590 - 378.004/A°-4 
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Figure captions 

Fig.l The energy curve obtained in axial RMF calculations as a function of the deformation 
parameter, (3, for the 52 Te and 54Xe isotopes. Calculated points are connected by dashed 
curves to guide the eye. Note that the curves for ^^^Te and ^^^Te are shifted downward 
by 0.02 MeV to distinguish them from the other curves. 

Fig. 2 The energy surface in the plane of deformation parameters, j3 and 7, calculated in 
the RMF theory with triaxial deformation for nuclei in the range of Z = 50 — 58 and 
= 64 — 72, arranged in the form of the nuclear chart. The energy difference between 
the contours is 1 MeV in total binding energy. The energy minimum is marked by the 
black region, in which the energy difference is less than 1 MeV. 

Fig. 3 The energy surfaces obtained in the triaxial RMF calculation without and with pairing 
are shown for ^^'^Te, ^^^Xe and ^^"^Ba. The energy difference between the contours is 1 
MeV in total binding energy. The energy minimum is marked by the black region. 

Fig. 4 The proton deformation parameter (3, obtained from the RMF calculation with axial 
deformation, is shown as a function of the neutron number. The experimental data are 



extracted from Ref. 45 . 



Fig. 5 The energy surfaces obtained in the triaxial RMF calculation with the NLl and TMA 
parameter sets are shown for ^^^Ba. The energy difference between the contours is 1 MeV 
in total binding energy. The energy minimum is marked by the black region. 

Fig. 6 The total binding energy calculated using the RMF theory with axial deformation is 
shown as a function of the neutron number. The experimental data are extracted from 



Ref. 46 . 



Fig. 7 The charge radius calculated using the RMF theory with axial deformation is shown as 



a function of the neutron number. The experimental data are extracted from Ref. ^\ . 
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